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ENVIRONMENTAL EFFECTS ON GASEOUS DISKS OF THE VIRGO SPIRAL GALAXIES 1 

HlROYUKI NAKANISHI 1,2 , NARIO KUNO 2 ' 3 , YOSHIAKI SOFUE 1 , NAOKO SATO 4 ' 5 , NAOMASA NAKAI 6 , YASUHIRO SHIOYA 7 ' 8 , 
TOMOKA TOSAKI 9 ' 2 , SACHIKO ONODERA 1 , KAZUO SORAI 4 , FUMI EGUSA 1 , AND AKIHIKO HlROTA 10 ' 2 

ABSTRACT 

We found high molecular fractions (/ mo i; ratio of the molecular to total gas surface densities) in three of five 
Virgo spiral galaxies in spite of their low total gas column density, based on 12 CO(/ =1—0) observations with 
the Nobeyama 45 m telescope equipped with a multi-beam receiver, BEARS. We interpret this as a result of 
environmental effects. Combining the CO data with HI data, the relationship between the surface density of the 
total gas (HI plus H2) and / mo ] indicates that the three galaxies near the cluster center have larger / mo i values 
than expected for field galaxies, while the others show normal / mo i. The large / mo i is interpreted as being due 
either to effective HI gas stripping, even in the inner disks, or to large ISM pressure induced by the high ICM 
pressure and/or ram pressure, although the possibility of an unusually high metallicity cannot be ruled out. 
Subject headings: galaxies: clusters: individual (Virgo), galaxies: ISM, ISM: molecules, radio lines: ISM 



1. INTRODUCTION 

Galaxy members in a cluster of galaxies are embedded in a 
hot diffuse intra-cluster medium (ICM) emitting X-rays, and 
suffering various environmental effects. Since Virgo is one of 
the nearest clusters, we can study the environmental effects 
there in greatest detail. 

Observation in the 21 cm line show that the HI gas is de- 
ficient for galaxies near the center of the Virgo Cluster (e.g., 
Cayatte et al. 1990). This is inter preted as be ing a result of 
ram pressure stripping. Vollmeretal. (2001) simulated the 
ram pressure effect by a numerical calculation and matched 
well t he observed distor t ion of the HI disk due to the ram pres- 
sure. Kennev & Young ( 1988, 1989) observed the Virgo spi- 
ral galaxies in the CO(7 =1-0) line using the 14 m telescope 
of the Five College Radio Astronomy Observatory (FCRAO) 
and found that the molecular gas content appears to be nor- 
mal, even in HI deficient galaxies. They reported that HI- 
deficient galaxies have a high ratio of CO flux to HI flux. 

Whether a galaxy suffers environmental effects or not, the 
inner disk of a spiral galaxy is dominated by molecular gas, 
because in general the inner disk is originally dominated by 
molecular gas. Therefore, it has not been clarified whether the 
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inner disk is affected by environmental effects or not. 

We can enumerate three possibilities for the physical condi- 
tion of the inner disk: (1) ram pressure stripping occurs only 
in the outer HI disk, and the inner disk suffers no environmen- 
tal effect, (2) ram pressure stripping occurs in the inner disk as 
well as the outer disk, but only HI gas is selectively stripped, 
(3) the ISM pressure of the inner disk increases due to exter- 
nal pressures, and molecular gas formation is enhanced. The 
molecular gas can be dominant in the inner disk in any of 
these three possibilities. 

In order to develop our understanding of the physical con- 
dition of the gas in inner disks of cluster spiral galaxies, we 
must pay attention to the molecular fraction (/moi; ratio of 
molecular to total gas surface density). In the first possibility, 
there would be no difference in the molecular fraction, / m oi> 
in the inner disk between cluster and field galaxies. However, 
the / mo i of cluster galaxies would be larger than that of field 
galaxies in the second and third possibilities. 

In this paper, we show results of the highest resolution 
single-dish observations in the CO line of five Virgo Clus- 
ter galaxies achieved with the Nobeyama 45 m telescope 
equipped with a multi-beam receiver, BEARS. We also com- 
pare the CO data with HI data with almost the same reso- 
lution obtained with the Very Large Array (VLA) C and D- 
configurations. These data sets enable us to investigate the 
molecular fraction at each point of the spiral galaxies with 
a fine scale, while earlier research has dealt with the total 
amounts of HI and CO (H2) gases in galaxies. Based on clues 
from /moi, we quantitatively discuss the physical condition of 
the gas disks in terms of (1) HI stripping, (2) external pressure 
due to the ICM or ram pressure, (3) intrinsic metallicity in a 
galaxy, and (4) the UV field of a galax y. The distance of th e 
Virgo cluster is taken to be 16.1 Mpc (Fere arese et all l996). 

2. OBSERVATIONS AND DATA 

The 12 CO(/ =1—0) observations of the Virgo Cluster spi- 
rals were made from 2002 December to 2004 April with the 
Nobeyama 45 m radio telescope. The half-power beam-width 
was 15" (1.2 kpc at the Virgo Cluster) at 1 15 GHz. We used a 
focal plane array receiver, BEARS (SIS 25-BEam Array Re- 
ceiver System), which consisted-of 25 beams (5 x 5) and was 
opera ted in the double-side-band (DSB) mode dSunada et alJ 
2000). The separation between the beams was 41 "1 . Twenty- 
five digital auto-correlators with a 500 kHz resolution (1.3 km 
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s~ J at 115 GHz) a nd a 512 MHz coverage (1332 km s" 1 ) 
( Sora i et all200 0a) were used as spectrometers. 

Calibration of the line intensity was made with an absorb- 
ing chopper wheel in front of the receiver, which yielded the 
antenna temperature T£ corrected for the atmospheric and an- 
tenna ohmic losses. Since the scale of the intensity varies 
with the beams, mainly due to the difference of the side band 
ratio between the upper side band (USB) and the lower side 
band (LSB), we scaled the observed DSB intensity to the SSB 
(single side band) intensity using scaling factors that were de- 
termined by comparing the CO intensities of NGC 7538 mea- 
sured with BEARS and a single-beam SSB receiver, SIS 100. 
The main beam efficiency of the telescope was T|mb = 0.40 
at 1 15 GHz, and the intensity given in this paper is the main 
beam brightness temperature, defined by Tmb = TaAImb- 

We mapped five Virgo galaxies with a grid spacing of 10"3. 
The number of total observed points for each galaxy was 576 
and the mapped area was 3 '.95 x 3 '.95, except for NGC 4254 
and NGC 4569, whose number of observed points was 1056 
and the area was 3'.95 x 7f38. The observations were made 
in the position-switching mode with an off-position at an off- 
set of 7' from the centers of the galaxies in azimuth. The 
typical rms noise of the main beam brightness temperature 
was ATmb = 0.09 — 0.18 K per velocity channel. Pointing of 
the antenna was calibrated by observing the continuum point 
source 3C273 at 43 GHz, and its typical error was less than 5" 
(peak valu e). Our observed sa mples were selected out of the 
galaxies in lSofue et al.l (12003). and the observing parameters 
are listed in Table 1. 

— Table 1 — 

To compare with the distribution of CO, we used the HI 
data of the target gala xies mapped with the VLA. The HI 
data were adopted from P hookun Vogel. & Mundvl ill 9931) for 
NGC 4254, from lCavatteetalJ (11990) for N GC 4402 NGC 
4569, and NGC 4579, and from lPhookun & Mundvl (119951) 
for NGC 4654. The spatial resolution of the HI maps was 
about 20", which is close to that of our CO maps. Table 2 
summarizes the parameters of the HI data. 

— Table 2 — 

Figure 1 shows the distribution of the CO integrated inten- 
sity, Iqo = J TwiBdv (K km s~'), (color) and the HI intensity 
(contours) of the five galaxies. The CO emissions distribute 
in the inner disks, and the extensions are typically up to radii 
of- l'.5 (7kpc). 

The HI gas in NGC 4254 and NGC 4654 is abundant and its 
distribution is strongly distorted. The HI disks of NGC 4402 
and NGC 4569 are truncated at the edges of the CO disks. In 
NGC 4579, the HI gas surrounds the central CO disk where 
the HI gas is deficient. 

— Figure 1 — 



3. ENVIRONMENTAL EFFECTS ON GASEOUS DISKS 

3.1. Molecular Fraction as a Function of the Surface Density 
of the Total Gas 

The fraction of molecular gas to the total gas (HI gas plus 
H2 gas), / mo i, is determined by the ISM (interstellar medium) 
pressure P, metallicity Z, and UV radiation U (Elmegreen 



1993). Here the molecular fraction / mo i is defined as 

Eh? 



/moi — 



-HI + Z -H 2 



(1) 



with Ehi and Eh 2 being the surface mass densities of HI 
and H2, respectively. Since the ISM pressure is approxi- 
mately proportional to the square of the gas surface density 
(Elmegree rIT993T) . 

P_ = E2 = (Ehi + Eh 2 ) 2 
Po Eg (E H i+Eh 2 )§' 

/ mo i i s expressed as a function of E, when Z and U ar e 
give n. iHonma. Sofue. & Arimotol ill 9951) . iKuno et al J (119951) . 
and ISorai et al. IHoOOb) investigated / mo i in nearby spiral 
gal axies, and showed that it is well reproduced by a model 
of Elmegreen (1993). Figure 2 shows model cu rves o f 
the E-/ mo i relation calculated based on Elmegr eenl (119931) . 
where we scale P, U, and Z with the values of the solar 
neighborhood: Po, Uq, and Zq. The gas surface density 
at the solar neighborhood, En, is t aken to be 8 M pc~ 2 
(Sanders, Solo mon. & Scovillell984l) . which is adopted to be 
1/2 

proportional to P . We show model curves taking UV radia- 
tions of 0.1£/o, lC/o, and 10[/n and metallicities of O.IZo, lZn, 
and lOZo. 

— Figure 2 — 

Figure 3 shows / mo i against E for each observed point 
of the five galaxies, where the CO-to-H2 conversion factor 



was adopted to be Xc r> = 1.0 x 10 /u H 2 cm~ z K -1 kirr 1 s 
( iNakai & Kunoll995l) . The dashed curves in all panels denote 
the averaged / mo i of NGC 4254 and NGC 4654. The molecu- 
lar fractions, f mo \, of NGC 4254 and NGC 4654 show similar 
curves to the models with U — IUq and Z = lZn, which is 
expected for field galaxies. 

On the contrary, the other three galaxies (NGC 4402, NGC 
4569, and NGC 4579) present larger molecular fractions than 
the dashed curve. NGC 4402 shows an extraordinarily large 
/moi in spite of the small surface density of the gas. For NGC 
4569 the plotted points are scattered on the upper side of the 
dashed curve. NGC 4579 shows a very large / mo i in spite of 
the small surface density of the gas as NGC 4402. In the three 
galaxies, f mo \ is always larger than the dashed curve. 

— Figure 3 — 

3.2. Origin of the Unusual Molecular Fraction in the Virgo 

Spirals 

3.2.1. Ram pressure Stripping 

Ram pressure stripping is a possible interpretation of the 
unusually large / mol of NGC 4402, NGC 4569, and NGC 
4579. 

Ram pressure stripping would occur if the criterion 

2 

> 1 (3) 



PlCM 



PlSM 



8v. 

V ro t 



is satisfied, where Pism and Picm are the volume number den- 
sity of the ISM and the ICM, respectively, R is the galactocen- 
tric radius of the element, d is the thickness of the gas disk, 8v 
is the relative velocity of an ISM cloud against the ICM, and 
V Iot is the rotational velocity of the galaxy. 

For a simple case, we consider a gas disk with thickness 
d — 0.05 kpc, rotating at V mt — 200 km s . For a typical 
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ICM wind with Pic m ~ 10~ 4 H cm' 1 and 5v ~ 1000 km s" 1 
in the Virgo Cluster (Hidaka & Sofue 2002), we obtain for the 
criterion of where ram pressure stripping can occur, 

p ISM < °- 05 (^) [Hem- 3 ]. (4) 

The HI gas can be dominant in the inter-arm region even 
in inner disks, and its density typically ranges from 0.01 to 1 
H cm -3 . Hence, the HI gas can be stripped even in the in- 
ner disks. On the other hand, the molecular cloud can hardly 
be stripped because the gas density in the molecular cloud is 
10— 1000 H2 cm -3 . Therefore, the HI gas can be stripped 
selectively from a disk, although the molecular gas still re- 
mains in the disk. In this case, f mo \ in the E — / mo i diagram 
moves toward the upper-left side of the E — f mo \ curve of field 
galaxies. 

NGC 4402, NGC 4569, and NGC 4579 are located within 
the projected radius of 2° from M87 at the center of the clus- 
ter, and are known to be HI deficient. Hence, they must have 
experienced ram pressure stripping of the HI gas. We empha- 
size that in these three galaxies, (1) the HI disks are restricted 
within the central CO disks of \ '.5 radii (7 kpc) (figure 1), (2) 
the HI clouds can satisfy the criterion of ram pressure strip- 
ping while the molecular clouds cannot, and (3) the E — f mo \ 
diagrams show unusually large / mo i, indicating that the in- 
ner disks are highly HI deficient (figure 3). These three facts 
support the scenario that the ram pressure stripping of HI gas 
occurred even in the inner disks. 

On the other hand, NGC 4254 and NGC 4654 are located 
as far as ~ 3° from the center of the cluster. They show the 
HI distribution extending to the outer disks and the E — / mo ] 
relation being the same as that of field galaxies. Hence, most 
of the HI gas must have not yet been stripped in spite of the 
strongly distorted HI disks. 

3.2.2. Higher External Pressure due to the ICM-Pressure or 
Ram Pressure 

Another possibility for making f mo \ large is a higher ex- 
ternal pressure due to the ICM-pressure, or the ram pressure, 
which might make the ISM pressure larger than that estimated 
by the surface density. 

The ICM pressure affects all gas disks isotropically. The 
typical density and temperature of the ICM at the Virgo 
Cluster are p ~ 10~ 4 cm~ 3 and T ~ 10 7 K, respectively 
llNulsen & B ohringerl ll995l) . On the other hand, the typical 
density and temperature of the interstellar molecular gas are 
p ~ 10 2 cm~ 3 and T ~ 10 K, respectively. Since the pressure 
P is proportional to pT, the ICM pressure and ISM pressure 
are comparable to each other near the cluster center. 

The ram pressure would also increase the ISM pressure if 
the ISM were not stripped, as mentioned by Kennev & Young 
( 1989). The typical velocity of a galaxy is ~ 1000 km s" 1 
near the Virgo Cluster center. The velocity dispersion of a 
molecular cloud is typically ~ 1 km s . Since the pressure 
P is proportional to pv 2 , the ram and the ISM pressures are 
comparable near the cluster center. As a result, the ISM pres- 
sure would increase and the higher ISM pressure would make 
/moi larger. 

We present the E — / mo i relation in figure 0] adopting 
(P/Pq) = (E/E0) 2 + n 2 . The second term indicates an incre- 
ment of the ISM pressure due to the ICM or the ram pressure. 
Figure|4]shows that the ICM pressure makes / mo i much larger. 
Thus, we can conclude that a higher external pressure due to 
ICM or ram pressures could be the origin of the large / mo i- 



3.2.3. Larger Metallicity 

Because the metallicity strongly affects / mo i, the larger 
metallicity might be an origin of the large / mo i. In order to 
examine whether the large metallicity makes f mo \ large, data 
on the metallicity distributions are necessary. However, there 
is no available data for these galaxies, except for NGC 4254. 
Therefore, we below present a new idea to calculate the metal- 
licity distribution using the CO, HI, and Ha data instead of 
using the metallicity data. 

The metallicity is expressed using the oxygen abundance 
(12 + log |Q/H] ) , which is correlated with the conversion fac- 
tor X ( Arimoto et all!996l) . We here define 

logZ=(12 + log[0/H]), (5) 

and adopt the relationship 

l0g (^) =_1 ° gZ + 93 - ^ 
The metallicity of the solar neighborhood logZo is taken to be 
8.9. We can calculate the H2 surface density using X if Z is 
given. Combining the HI and H2 surface densities, we obtain 
the total gas density E, which gives the pressure P based on 
equation (|2jl. 

The UV strength can be calculated using the Ha data. Ha 
images of the Vi rgo Cluster galaxies, except for NGC 4402, 
were archived bv Koonm ann. Kennev. & Young ll2001). Fig- 
ure [5] shows contours of Ha superimposed onto the CO im- 
ages. We convolved the Ha image with the same beams as 
the HI data. We adopted U = 2.42 x 10~ 7 erg cirr 2 s -1 sr _1 
jRevnoldsll984l) . For NGC 4402, we estimated U while con- 
sidering the Schmidt law, 

where we adopted n = 1.33 (Komu giet alJ20 05). 

Thus, we can calculate X, P, and U if Z is given. Us- 
ing obtained P, U, and Z, we can calculat e the mol ecular 
fraction / mo i based on the model of Elmegr eenl ( 119931) . On 
the other hand, we can independently calculate f mo \ using 
only the HI and H2 surface densities. We search for an ap- 
propriate metallicity, which gives the same / mo i in the two 
ways. The searching range of the metallicity is adopted to be 
6.68 <logZ< 11.12. 
This procedure is presented as a flow-chart in figure|5] 
We plotted the calculated metallicity against the radius for 
each galaxy (figure 0. Each point corresponds to each ob- 
served point. The metallicities of NGC 4254 and NGC 4654 
show that logZ is about 9.5 - 10 at the galactic center, and that 
it graduall y declines. The metallici t y of N GC 4254 was cal- 
culated bv lVila-Costas & Edmunds! i 19921) and Zaritsk v" et alJ 
( 119941) . Our obtained values are almost consistent with their 
results. 

On the other hand, the metallicities logZ of NGC 4402, 
NGC 4569, and NGC 4579 often exceed 10. However, 
the former research shows that there is no galaxy that gives 
logZ > 10. In addition, the dispersions of logZ of NGC 
4402, NGC 4569, and NGC 4579 are larger than the oth- 
ers. Moreover, the metallicities logZ of NGC 4402 and NGC 
4569 increase with the radius. This tendency is unnatural be- 
cause t he metal licity logZ usually decreases with th e radius 
dVila-Costas & Edmundsll 19921 IZaritskv et ai1 ll994l) . In the 
case of NGC 4579, there are many points where Z cannot be 
calculated within the solution range 6.68 < logZ < 11.12 in 
the inner region. 
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Table 1 
Observational Parameter 



NGC R.A.(J2000) Dec.(J2000) Morph. B T i P.A. V hcl 

h m s d m s mag. deg. deg. kms~' 

4254 12 18 49.61 +14 24 59.6 SA(s)c 10.44 42 68 2405 

4402 12 26 07.45 +13 06 44.7 Sb 12.55 75 90 234 

4569 12 36 49.82 +13 09 45.8 SAB(rs)ab 10.26 64 23 -235 

4579 12 37 43.53 +11 49 05.5 SAB(rs)b 10.48 38 05 1520 

4654 12 43 56.67 +13 07 36.1 SAB(rs)cd 11.10 52 125 1039 

Col. (1): Galaxy name. Cols. (2) and (3): Central position taken from Sofue et al. 1 2003). Col. (4): Morphological type from RC3 i de Vaucouleurs et al. 1991). 

Col. (5): Tota l B-band magnitude taken from RC 3. Cols. (6) and ( 7): The inclination angle and the position angle (P.A.) taken from 

Koopmann, Kennev. & Youna 1 2001 ) or Phookun. Voael, & Mundv 1 1993). Col. (8): Systemic velocity taken from Kennev & Youna 1 1988) or RC3. 



Therefore, the possibility of higher metallicity is less plau- 
sible, because an abnormal metallicity distribution is neces- 
sary to reproduce such a large f mo \. However, this possibility 
still cannot be ruled out, since nobody has directly measured 
the metallicities of NGC 4402, NGC 4569, and NGC 4579. 

3.2.4. Lower UV Radiation Field 

Because the UV strength affects / mo i, a lower UV radiation 
field might be the origin of the large f mo \. In order to examine 
this effect, we plotted the relation between the UV strength 
and the total gas density (figure[S}, excluding the central 1 kpc 
to avoid Ha emission from active galactic nuclei. This figure 
shows that there is little difference in the UV strength, U /Uq, 
for E < 20M Q pc~ 2 , although U/U is lower in NGC 4569 
than in NGC 4254 and NGC 4654 in the case of E > 2OM 
pc~ 2 . Thus, it is difficult to attribute the large f mo \ to a lower 
UV radiation field, because a large / mo i for E < 20M© pc -2 
cannot be explained by a difference in the UV strength. 



4. SUMMARY 

We observed five Virgo spiral galaxies with the NRO 45 m 
telescope with BEARS in the 12 CO(7 =1-0) line. Com- 
paring the CO data with the HI data to investigate the envi- 
ronmental effect, we found that (1) the HI gas disks of NGC 
4402, NGC 4569, and NGC 4579, located near the center of 
the cluster, are restricted within the inner disks, (2) the HI 
clouds can satisfy the criterion of the ram pressure stripping, 
while the molecular clouds cannot, and (3) the E-/ mo i dia- 
grams of the three galaxies show an unusually large / mo i for 



a low total gas column density. This large f mo \ might imply 
that the HI gas is selectively stripped. On the other hand, NGC 
4254 and NGC 4654, located far from the cluster center, have 
extended HI disks and show the normal / mo i expected for field 
galaxies. Most of the HI gas of these two galaxies must not 
have been stripped, in spite of the strongly distorted HI disks. 

A higher external pressure due to the ICM-pressure and/or 
ram pressure might be another possibility causing the unusu- 
ally large f moh 

Moreover, we examined whether differences in the intrin- 
sic conditions (metallicity and UV) cause the unusually large 
/ mo l. As a result, we found that an unusually high metallic- 
ity is necessary to explain such a large / mo i, although it still 
cannot be ruled out. 

We also found that the large f mo \ cannot be explained by 
only a difference in the UV radiation field. 

From these discussions, we note that the E — / mo ] diagram 
is a good tool for investigating the environmental effect. An 
abnormal f mo \ indicates that the ISM suffers cluster environ- 
mental effects: ram pressure stripping or higher external pres- 
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Table 2 
HI Data Parameter 



NGC Synth. Beam Size AV AT B Reference 





arcsec 2 


km s 


K 




4254 


25" ) 


< 24" 


10.3 


0.45 


1 


4402 


21" ) 


< 17" 


25.0 


1.7 


2 


4569 


17" > 


< 13" 


25.0 


2.8 


2 


4579 


19" ; 


< 18" 


25.0 


2.7 


2 


4654 


24" ) 


< 25" 


10.3 


0.44 


3 



Col.(l): Galaxy name. Col. (2): Synthesized beam size. Col. (3): Velocity resolution. Col. (4): R.m.s. noise in brightness temperature. Col. (5) 
References.— (l lPhookun. Vogel. & Mund^fi993T) : (2 1Cavatte et"aTHT^Cl) : G IPhookun & Mund\l ft995) 
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FIG. 1 . — Integrated CO intensity maps (color) and the integrated HI intensity maps (contour). The contour levels are 100, 200, 400, 800, 1600 K km s 
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FIG. 3. — Molecular fraction, f mo \, against the surface density of the total gas E. The dashed curves in all panels denote the averaged / mo i of NGC 4254 and 
NGC 4654. The data of E K > 3.8 x 10 20 H cm" 2 or Eh 2 > 7.0 x 10 19 H 2 cm" 2 (3a) are plotted. 
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